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Current-driven  electrostatic  ion-cyclotron  instability  has  so  far  been  studied  for  Maxwellian  plasma 
with  isotropic  and  anisotropic  temperatures.  Since  satellite-measured  particle  velocity  distributions 
in  space  are  often  better  modeled  by  the  generalized  Lorentzian  (kappa)  distributions  and  since 
temperature  anisotropy  is  quite  common  in  space  plasmas,  theoretical  analysis  of  the  current- 
driven,  electrostatic  ion-cyclotron  instability  is  carried  out  in  this  paper  for  electron-proton  plasma 
with  anisotropic  temperatures,  where  the  particle  parallel  velocity  distributions  are  modeled  by 
kappa  distributions  and  the  perpendicular  velocity  distributions  are  modeled  by  Maxwellian  distri¬ 
butions.  Stability  properties  of  the  excited  ion  cyclotron  modes  and,  in  particular,  their  dependence 
on  electron  to  ion  temperature  ratio  and  ion  temperature  anisotropy  are  presented  in  more  detail. 

For  comparison,  the  corresponding  results  for  bi-Maxwellian  plasma  are  also  presented.  Although 
the  stability  properties  of  the  ion  cyclotron  modes  in  the  two  types  of  plasmas  are  qualitatively  sim¬ 
ilar,  significant  quantitative  differences  can  arise  depending  on  the  values  of  Ke  and  k,-.  The  com¬ 
parative  study  is  based  on  the  numerical  solutions  of  the  respective  linear  dispersion  relations. 

Quasilinear  estimates  of  the  resonant  ion  heating  rates  due  to  ion-cyclotron  turbulence  in  the  two 
types  of  plasma  are  also  presented  for  comparison.  ©  2011  American  Institute  of  Physics. 

[doi:  10. 1063/1 .3632974] 


I.  INTRODUCTION 

The  satellite-measured  particle  velocity  distributions  in 
the  solar  wind  and  in  many  space  plasmas  often  exhibit  non- 
Maxwellian  suprathermal  tails  that  decrease  as  a  power-law 
of  the  velocity. 1  Such  deviations  from  the  Maxwellian  distri¬ 
butions  are  expected  in  plasmas  with  sufficiently  low  degree 
of  collisionality.  The  distribution  function  that  can  better 
model  such  particle  velocity  distributions  is  the  so-called 
generalized  Lorentzian  or  kappa  (k)  distribution.2  The  kappa 
distribution  with  a  finite  value  of  the  spectral  index  k  has  a 
power-law  tail  at  velocities  larger  than  the  thermal  velocity 
and,  consequently,  has  a  substantially  larger  number  of 
suprathermal  particles  in  comparison  with  the  Maxwellian 
distribution.  It  approaches  the  Maxwellian  distribution  in  the 
limit  as  k  — >  oc.  Typical  values  of  k  for  space  plasmas  are  in 
the  range  2-6.  In  the  last  several  years,  many  authors  studied 
electrostatic  and  electromagnetic  waves  in  spatially  homoge¬ 
neous  and  weakly  inhomogeneous  magnetoplasmas  using 
different  types  of  kappa  distributions  for  the  equilibrium 
state. 1,3-17  The  presence  of  substantially  larger  number  of 
suprathermal  particles  in  comparison  with  the  Maxwellian 
was  shown  to  have  important  distinguishing  effects  on  the 
spectral  properties  of  the  excited  waves  and  on  the  wave- 
particle  interactions,  in  particular.  In  this  paper,  we  present 
for  the  first  time  the  stability  properties  of  the  electrostatic 
ion  cyclotron  modes  in  current-carrying  Lorentzian  (kappa) 
plasma  with  anisotropic  electron  and  ion  temperatures  and 
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compare  them  with  those  in  Maxwellian  plasma.  The  com¬ 
parative  study  is  based  on  the  numerical  solutions  of  the  lin¬ 
ear  dispersion  relations  for  the  two  types  of  plasma. 
Numerical  analysis  allows  the  investigation  of  parameter 
regimes  that  are  not  accessible  to  tractable  analytical  treat¬ 
ments  and  thus  allows  a  more  comprehensive  quantitative 
comparison.  Quasilinear  estimates  of  the  resonant  ion  heat¬ 
ing  rates  due  to  ion  cyclotron  turbulence  in  the  two  types  of 
plasma  are  also  presented  for  comparison. 

Electrostatic  ion-cyclotron  instability,  which  considers 
waves  propagating  at  large  angles  to  the  ambient  magnetic 
field  with  frequency  near  the  ion  cyclotron  frequency,  can 
arise  from  various  free  energy  sources  such  as  field-aligned 
currents,  ~  electron  beams,  ion  beams  ’  or  combined 
effects  of  ion  beams,  and  counterstreaming  electrons.27  We 
concentrate  here  on  the  field-aligned  current,  which  is  a  com¬ 
mon  feature  in  many  space  and  laboratory  plasmas,  as  the 
free  energy  source  for  the  ion-cyclotron  instability. 

Current-driven  electrostatic  ion-cyclotron  instability 
plays  an  important  role  in  the  generation  of  ion  cyclotron  tur¬ 
bulence  and  in  the  concomitant  anomalous  ion  heating 
observed  in  space  and  laboratory  plasmas.  The  previous  the¬ 
oretical  and  numerical  studies  of  the  instability  assumed  that 
the  electron  and  the  ion  distributions  in  velocity  space  are 
Maxwellians.  Drummond  and  Rosenbluth 18  first  developed 
the  theory  of  this  instability  for  Maxwellian  plasma  with  iso¬ 
tropic  electron  and  ion  temperatures  and  showed  that,  for 
Te  rs  Tj,  the  critical  value  of  the  electron  drift  relative  to  the 
stationary  ions  (i.e.,  the  critical  current)  for  the  onset  of  the 
instability  is  much  smaller  than  that  for  the  onset  of  the  elec¬ 
trostatic  ion-acoustic  instability,  which  considers  waves 
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propagating  parallel  to  the  magnetic  field.  Kindel  and  Ken¬ 
nel19  extended  the  theoretical  analysis  of  Drummond  and 
Rosenbluth 1  ’  and  showed  that  the  ion  cyclotron  waves  are 
unstable  to  smaller  currents  for  a  broad  range  of  Te/Tj.  The 
effects  of  temperature  anisotropies,  which  are  also  important 
features  of  collisionless  plasmas,  on  the  excitation  of  the 
current-driven  electrostatic  ion-cyclotron  instability  in  Max¬ 
wellian  plasma  were  first  investigated  by  Lee,21  who  fol¬ 
lowed  the  theoretical  analysis  of  Kindel  and  Kennel.19  Later, 
Okuda  and  Ashour-Abdalla22'23  investigated  the  effects  of 
temperature  anisotropies  over  a  wider  range  of  parameters 
by  numerically  solving  the  dispersion  relation.  For  a  review 
of  these  previous  studies,  see  Ref.  28. 

In  Sec.  II,  we  present  the  mathematical  model  and  the  lin¬ 
ear  dispersion  relations  for  the  current-driven  ion  cyclotron 
modes  in  Maxwellian  and  Lorentzian  (kappa)  plasmas.  In 
Sec.  Ill,  we  present  the  approximate  analysis  of  the  dispersion 
relations  in  order  to  gain  some  preliminary  understanding  of 
the  differences  in  the  stability  properties  of  the  ion  cyclotron 
modes  in  the  two  types  of  plasma.  In  Sec.  IV,  we  present  the 
numerical  solutions  of  the  dispersion  relations.  In  Sec.  V,  we 
calculate  the  resonant  ion  heating  rates  due  to  ion-cyclotron 
turbulence  in  the  two  types  of  plasma,  within  the  framework 
of  the  quasilinear  theory,  and  discuss  their  differences.  The 
paper  is  concluded  with  a  summary  in  Sec.  VI. 


for  I m  at  >  0.  Here  Fa(v2 ,  v  j  is  the  unperturbed  (equilib¬ 
rium)  distribution  function  for  the  charged  particle  species 
a  =  e,i,  Q  =  eBo/(mjc),  Jn(fF)  is  the  Bessel  function  of 
order  n,  /(  =  kxvx/Q and  k2  =  k2  +  k\.  For  analytic  con¬ 
tinuation  to  Im  co  <  0,  the  Landau  contour  has  to  be  used  for 
carrying  out  the  vy  integration.  The  two  types  of  equilibrium 
distribution  function,  considered  here  for  the  derivation  of 
the  linear  dispersion  relation,  are  the  following. 

A.  Bi-Maxwellian 
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where  uay  =  v  -  Vao  with  V,/(l  =  Vo  for  the  electrons  and 
Vco  =  0  for  the  ions,  92 1  (  is  related  to  the  particle  tempera¬ 
ture  raj|(x)  by  02y(x)  =  27;|| (x)/»ia  =  2V2||(1),  and  Fx  is  nor¬ 
malized  to  the  particle  density  «o(same  for  both  electrons 
and  ions).  The  definitions  of  Ta y  and  Ta  are 


n0Ta  ||  =  2nmx 
noTx±  =  nmc 


<7vX(7v||VXu2||F’0!(v2  ,  V||),  (3) 

fifvX(7v||  v3±Fx  (v2  ,  V|| ) .  (4) 


II.  LINEAR  DISPERSION  RELATION  FOR 
CURRENT-DRIVEN  ION  CYCLOTRON  MODES 

The  starting  point  of  our  study  is  the  linear  dispersion 
relation  for  obliquely  propagating  electrostatic  ion  cyclotron 
modes  in  current-carrying  plasma  with  anisotropic  electron 
and  ion  temperatures.  For  this,  we  consider  spatially  homo¬ 
geneous,  nonrelativistic,  collisionless  plasma  in  which  elec¬ 
trons  are  drifting  along  the  uniform  ambient  magnetic  field 
Bo  with  velocity  Vo  relative  to  the  stationary  ions.  Consider¬ 
ing  normal  modes  of  the  form  exp  ;(k  ■  r  —  cot)]  with 
frequency  co  and  propagation  vector  k,  we  assume 
| (co  -  k||Vo)/Qt.|  <C  L  |fc||Ve||/i2f|  1  and  k±Ve±  /Qe  <C  1, 

where  Qe  =  eBo/(mec)  is  the  electron  cyclotron  frequency, 
Vf||  and  Ve±  are  the  electron  thermal  speeds  defined  as 
V2|,(±)  =  Ten (X)/»ie,  kj|  =  (k  •  Bo)/fio.  and  k±  is  the  perpen¬ 
dicular  (to  Bo)  component  of  k.  We  further  assume  that 
cope/Qe  <C  1,  where  cope  is  the  electron  plasma  frequency. 
These  assumptions  imply  that  the  electrons  are  strongly  mag¬ 
netized  particles  so  that  their  dynamic  response  to  the  elec¬ 
tric  field  perturbation  is  practically  one-dimensional  (parallel 
to  Bo).  For  the  ion  response,  however,  we  retain  both  the  par¬ 
allel  and  the  perpendicular  dynamics.  Then,  according  to  the 
linearized  Vlasov  theory,  the  dispersion  relation  for 
obliquely  propagating  electrostatic  waves  is  given  by29 


Substituting  /^(v2  ,V|j),  given  by  Eq.  (2),  into  Eq.  (1)  and 
performing  the  velocity  space  integrations,  we  obtain  the 
well-known  dispersion  relation 
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Here,  co2x  =  4ne2n0/mx,  A „(&,•)  =/„(&,-)  exp(-6,-),  /„  is  the 
modified  Bessel  function  of  the  first  kind,  b,  =  kx02x/  (2Q2), 
Qe  =  (co-  k||Vo)/(k||0(,||),  and  Qni  =  (co  -  nQ;)/(fc||0/||), 
where  k«  >  0  is  assumed.  Zm(q)  is  the  well-known  plasma 
dispersion  function^  associated  with  the  Maxwellian  vy -dis¬ 
tribution  and  it  is  defined  as 


Zm(q)  = 


■+0°  ,  exp(— s2) 
as - 

s-c 


(6) 


for  Im  q  >  0  and  as  the  analytic  continuation  of  this  for 
Img  <  0.  The  prime  notation  on  denotes  its  derivative 

with  respect  to  the  argument. 


B.  Kappa-Maxwellian 
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where /(kb)  =  r(fca  +  l)/[/cy2r(?ca  +  1/2)],  T(x)  being  the 
gamma  function.  As  before,  uay  =  vy  —  V^o  with  Vao  =  Vo 
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for  the  electrons,  V7q  =  0  for  the  ions,  and  F 7  is  normalized 
to  the  particle  density  «o(same  for  both  electrons  and  ions), 
while  0^||  and  O2  are  related  to  the  particle  “temperatures” 
T4  and  Taj_,  as  defined  by  Eqs.  (3)  and  (4),  according 
to  6\ ||  =  [(2kx  -  l)/Ka\(T^/ma)  =  [(2 kx  -  1)/k«]V^|1  and 
d2±  =  2T7  i /my  =  2U2±,  provided  Ky  >  1/2.  As  /c,  — >  oo, 
Fa  [given  by  Eq.  (7)]  asymptotically  approaches  the  bi- 
Maxwellian  [given  by  Eq.  (2)].  In  Eq.  (7),  the  particle  distri¬ 
bution  in  V|| — space  is  modeled  by  the  kappa  distribution 
function,  while  the  distribution  in  v  — space  is  modeled  by 
the  Maxwellian  distribution.  The  physical  arguments  for  not 
considering  a  kappa  distribution  in  v  — space  are  the  follow¬ 
ing.  First,  in  the  presence  of  an  ambient  magnetic  field,  the 
mechanism  that  produces  the  kappa  type  distribution  in  ve¬ 
locity  space  is  most  likely  to  be  more  effective  in  vy — space. 
Second,  the  most  important  physical  aspect  of  the  kappa  dis¬ 
tribution,  which  distinguishes  it  from  the  Maxwellian  distri¬ 
bution,  is  the  resonant  wave-particle  interaction  (Landau  and 
cyclotron  resonances)  between  the  excited  waves  and  the 
enhanced  population  of  suprathermal  charged  particles  that 
are  moving  along  the  ambient  magnetic  field.  It  may  be  men¬ 
tioned,  however,  that  a  kappa  distribution  in  v  — space  does 
not  pose  a  serious  mathematical  difficulty.  With  it,  the  veloc¬ 
ity  integral  in  the  ion  term  involving  J2  and  ( OF,  / <9v ,  )  can¬ 
not  be  expressed  in  terms  of  any  known  mathematical 
function;  but,  it  can  be  numerically  evaluated  quite  easily. 
The  integral  approaches  its  Maxwellian  counterpart,  A „(hj), 
in  the  limit  as  k ,•  — »  oo,  and  its  numerical  values  for  finite 
values  of  k,  are  not  too  different  from  those  of  A „{bi). 

Substituting  Fa(v2,  vn),  given  by  Eq.  (7),  into  Eq.  (1) 
and  performing  the  velocity  space  integrations,  the  disper¬ 
sion  relation  is  obtained  as 


III.  APPROXIMATE  STABILITY  ANALYSIS 

Here,  we  present  the  approximate  analytical  solutions  of 
the  dispersion  relations  in  order  to  gain  some  preliminary 
understanding  of  the  ion-cyclotron  instability  in  the  two 
types  of  plasmas.  Considering  co  =  Re  co  +  i  Im  co  and  mak¬ 
ing  a  Taylor  series  expansion  of  D(k,co)  around  co  =  Re  co 
while  assuming  Im  co  <C  Re  co  and  bn  D  <C  Re  /J,  we  find  to 
lowest  order, 

Re  D(k,  co  =  Re  co)  =  0,  (10) 

Im  D( k,  co  =  Re  co) 

lmco=-  .  1  ;  .  (11) 

[<9Re  D(k,  co)/dco\(0=Rem 

These  two  equations  determine  Re  co  and  Im  co.  Generally 
speaking,  unstable  solutions  are  found  when  Re  (o/k\\  <  Vo 
(i.e„  Re  <  0)  so  that  (dFe/dv||)V||=Reaj/jt||  >  0  (Landau 
growth)  and,  simultaneously,  when  (a)  Regm,  >  1  for  cyclo¬ 
tron  modes  with  Re  co  «  mQ,;  and  (b)  Reg„,  1  for  all 
n  ^  m,  such  that  the  Landau  growth  rate  due  to  the  drifting 
(current-carrying)  electrons  exceeds  the  total  ion  damping 
(both  Landau  and  cyclotron  damping)  rate.  For  the  approxi¬ 
mate  stability  analysis,  we  follow  the  commonly  used  analyt¬ 
ical  approach  of  Drummond  and  Rosenbluth,18  i.e.,  we 
consider  |;e|  <C  1  and  |  gni  |  1  for  all  n  and  retain  the  lead¬ 

ing  terms  in  the  power  series  expansion  of  the  electron 
plasma  dispersion  function  and  the  asymptotic  expansion  of 
the  ion  plasma  dispersion  function.  The  power  series  and  the 
asymptotic  expansions  of  ZM{g)  are’0 

ZM(g)  =  /V^exp(-g2)  -  2g  +  4g3/3  -  •  •  • ,  |g|  <  1,  (12) 


co2  co2-  +0° 

D(k,co)  =  l--^-Z'  (ge)~ 


k2o : 


k202 

K  ui ||  n=— oo 


z!  (c  )  -  2n^‘  0,11  z  (c  ■) 
k\\6i\\  Oll 


=  0, 


(8) 


ZM(?)  ^ /0texp(-g2)  |?|  >  1,  (13) 

?  2gJ 

while  the  power  series  and  the  asymptotic  expansions  of 
ZKa( s)  for  integer  values  of  Ky  are3'17'18 


where  A n(bi)  =  I„(bj)  exp(-bj),  I„  is  the  modified  Bessel 
function  of  the  first  kind,  b,  =  k\  02L/(2Fl2), 
Qe  =  (co-  &||Vo)/(£||0e||),  and  gni  =  (co  -  nQ)  /  (k^d^)  with 
^11  >  0.  The  dispersion  relation  for  kappa-Maxwellian 
plasma  is  formally  equivalent  to  that  for  the  bi-Maxwellian 
plasma,  except  that  02^  =  {Ty^/mr,)[{2Ky—  1)/k2]  and 
ZKa(g)  is  the  modified  plasma  dispersion  function  associated 
with  the  kappa  distribution  function  in  V|| — space.  The  modi¬ 
fied  plasma  dispersion  function  is  defined  by16'17 


ZKa(c)  =  iy/n 


f(Ka) 


(1  +g2/K^ 


2k„ 


1 


2k„ 


3  Kr/ 


Id  <  1, 

(14) 


ZKa{g)  =  iy/n 


f{Kg) 

(1  +  g2/Ky)K‘+1 


1  Kg  1 

?  2k«  -  1  ?3 


ld»L 

(15) 


2Ka(?)  = 


'  +  OO 


ds 


-00  (s-g)(l  +s2/kx) 


K«  +  1 


(9) 


for  Im  g  >  0  and  by  the  analytic  continuation  of  this  for 
Im  g  <  0.  The  prime  notation  on  ZKa  (5)  denotes  its  derivative 
with  respect  to  the  argument.  Aside  from  the  multiplication 
factor,  this  function  is  analogous  to  the  plasma  dispersion 
function  that  was  first  introduced  and  discussed  by  Summers 
and  Thorne3  and  was  later  further  analyzed  by  Mace  and 
Hellberg31  and  Summers  et  al. 32 


For  noninteger  (including  half-integer)  values  of  Ka,  the  as¬ 
ymptotic  expansion  has  small  correction  terms,31’32  which 
may  be  neglected. 

A.  Bi-Maxwellian  plasma 

Restricting  ourselves  to  the  fundamental  ion  cyclotron 
mode,  i.e.,  Re  co  «  Q,,  we  retain  only  the  n  =  I  term  in  the 
ion  sum.  Using  the  leading  terms  of  Eqs.  (12)  and  (13)  in 
Eq.  (5),  we  find 
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Re  Z)(k,  m  =  Re  m)  =  1 


He  Hd\  I 


a 


k2e24  k2e\ 


Ai  (bi), 
(16) 


,  ,  _  or  Rem  —  A'iiVo 

Im77(k,m  =  Rem)  =  2iJh — - — ^ — 

v  ’  y  V  ^ 2^2  7-  73 


HI 


^11 


+  2i>/7t 


ffl2  Rem-Q-^JQ, 


A20; 


*11 0/11 


x  A1(ft/)exp(-^,.). 


(17) 


Then,  referring  to  Eqs.  (10)  and  (11),  we  obtain 

Rem  _  Te\\  7^  A t(fe,-) 

~  7)||7)xl+A24e’ 


(18) 


1m  m 

”oT 


/nr,||  Till  Ai(fc;)  r^||  V0  —  Re  m 

V  2  7^721(1+^2  [  k[lv4 


Te\\  Rem  -  (1  -  Tq/Ti±)Qi 

T)\\  ^Vt|i 


A1(6i)exp(-?/?) 


(19) 


where  17,-  =  (Rem  —  Q,)/(v/2^|| V/||).  We  have  related  dl\\(±) 
to  7+(x)  and  introduced  the  Debye  length,  lDe,  defined  by 
)?De  =  T4/{me(o2e).  The  onset  condition  for  instability 
(Imm  >  0)  is 


Vo 

^11 


> 


Re  m 

We\\ 


+ 


Te\\  Rem  -  (1  -  r,1|/T/±)Q, 

^11  WH 


A1(h,)exp(-^). 


(20) 


These  results  for  anisotropic  ion  temperatures  (Tp  ^  T,  \ ) 
were  previously  obtained  by  Lee21  and  by  Okuda  and 
Ashour-Abdalla.2-  They  reduce  to  the  earlier  results  of 
Drummond  and  Rosenbluth18  for  isotropic  ion  temperatures 

(7)1,  =  TiL). 


B.  Kappa-Maxwellian  plasma 

As  before,  we  restrict  ourselves  to  Rem  «  Q,  and  use 
the  leading  terms  of  Eqs.  (14)  and  (15)  in  Eq.  (8)  to  find 


ReD(k,  m  =  Re  m)  =  1 


( tope  2  Kc  +  1 

w 


Hi  k 


Q, 


k2d2i\ |  92±  Re  co  -  Q, 


Ai(*,).  (2D 


ImZ)(k,m  =  Rem)  =2i\/n 


orpe  Rem  — A||V0kc  +  1 


*||  94 


f{.Ke) 


+  2iy/n 


03 Pi  Ai  (bi)  /(k,-) 

k%  kfy  (l+?2/Kl.)Ki+1 

IRem-Q,  ,  0|  X 


Kj  i  +  q,./fc,-  el 


(22) 


,A 


Then,  referring  to  Eqs.  (10)  and  (11),  we  obtain 

Re  m 
Q,  ~ 

Te II  ?)||  2?cg  -  1 _ Ai  (fe/) _ 

+  r,-||  Tj±  2 k„  +11  +  [(2 Ke  -  1)/(2jv  +  1  )]k2k2De  ’ 

(23) 


Imm 

Q, 


x 


x 


fn/2Ke-  iyrf||  T,-|| _ Ai(fe/) _ 

V2\2k,  +  l)  r,-||  7k  {1  +  [(2fCe  -  l)/(2ice+  l)]k2A2De}2 


A||V0  -  Remx-f  +  1 
A||V„||  Ke 


3/2 

f(*e) 


Te II  Ai(fc,-)  /  2 Ki  \ 
Ti\\  A||V,-||  \2iq-lJ 


3/2 

/(Ki) 


Vf 


Ki  -  1  /2 


-(«.'+!) 


K/-  1/2 


Ki  +  1 

Ki 


(Re  m 


Q/)  + 


2 ^  -  1  r,-n 

2k,-  7k 


(24) 


where  17,-  =  (Rem  —  fl,-)/(\/2A||V,j|).  We  have  related  ^II(-L) 
to  7km  and  introduced  the  Debye  length  XDe.  The  onset 
condition  for  instability  (Im  m  >  0)  is 


V0  Rem  T4 

Vel\>W^\\+fi\\ 

2k,  \ 3/2  / 2Ke-  1\ 3/2  Ke  f{Ki)Aj(bj) 
2k,  —  1/  \  2 Ke  J  Ke+lf(Ke)  k\\Vi\\ 

(«.-+!) 


X  1  + 


>17 


Ki  — 1/2 


1+ 


>17 


Ki-  1/2 


K,+  l.  2k,—  1  7’,-m 

x - (Re  m  —  Q, )  H — - ^-Q,- 


2k,-  Tl 


,'X 


(25) 


The  results  presented  in  Eqs.  (21 )— (25)  are  new  and  they 
reduce  to  the  corresponding  results  for  bi-Maxwellian 
plasma  [Eqs.  ( 1 6)— (20)]  when  Ke:  k,  — >  00. 

The  approximate  analytical  solutions  indicate  that  the 
stability  properties  of  the  current-driven  ion  cyclotron  modes 
in  both  types  of  plasma  (bi-Maxwellian  and  kappa- 
Maxwellian)  depend  on  me/mi,  (ope/kle,  Vo,  k\\,  k±,  electron 
to  ion  temperature  ratio  (Te«/T4),  and  ion  temperature  ani¬ 
sotropy  (T,-||/T, ■_)_).  Dependence  on  cope/(7le  arises  through 
the  term  a>pe / (A2  02ii)  in  the  dispersion  relations  as 
®pe/(^20e||)  (£24J_1  =  KM)(m^/Q2)/(A2p2),  where 

Pi(=  Vi±/Qj)  is  the  ion  gyroradius.  The  analytical  solutions 
also  indicate  that  Rem,  Imm,  and  the  onset  condition  for 
instability  in  kappa-Maxwellian  plasma  can  be  quite 
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FIG.  1 .  Threshold  (critical)  value  of  electron  drift  speed,  normalized  to  V), 
and  denoted  by  Vc,  as  a  function  of  rf||/r, ||  and  7)  |  /7)  i  for  a>pefCle  =  1/15. 
The  solid  curves  represent  results  for  bi-Maxwellian  plasma  and  the 
dashed  curves  represent  results  for  kappa-Maxwellian  (k6  =  k,  —  3) 
plasma.  The  curves  are  labeled  by  the  selected  values  of  T^/Tj±  and 

r«ll /in¬ 


different  from  those  in  bi-Maxwellian  plasma  depending  on 
the  choice  of  the  values  of  Ke  and  k,  .  Comparison  shows  that 
Re  cd  is  reduced  in  kappa-Maxwellian  plasma,  while  both  the 
electron  term  that  drives  the  instability  when  Vo  >  Re  co/k^ 
and  the  ion  cyclotron  damping  term  in  Im  cd  are  enhanced 
in  kappa-Maxwellian  plasma.  The  enhancement  of  the 
electron-drive  term  in  kappa-Maxwellian  plasma  is  due  to 
the  reduction  of  Re  cd  and  the  multiplying  factors  involving 
Ke.  The  enhancement  of  the  ion  cyclotron  damping  term  in 
kappa-Maxwellian  plasma  is  predominantly  due  to  its 
power-law  dependence  on  >  1  in  contrast  with  its  expo¬ 
nential  dependence  on  >  1  in  bi-Maxwellian  plasma.  For 
the  same  reason,  the  threshold  value  of  Vo  for  the  onset  of 
the  instability  in  kappa-Maxwellian  plasma  is  enhanced  over 
the  corresponding  value  in  bi-Maxwellian  plasma. 

We  shall  not  pursue  the  analytical  solutions  any  further 
as  they  are  valid  only  for  the  restricted  values  of  the  afore¬ 
mentioned  parameters  for  which  the  assumed  conditions 
(|?e,|  <C  1  and  c„,  2>  1  for  all  n)  can  be  satisfied.  Outside  the 
range  of  validity,  the  analytical  solutions  are  not  only  quanti¬ 
tatively  inaccurate;  they  even  give  erroneous  dependence  on 
k|| ,  for  example.  For  an  accurate  and  comprehensive  study  of 
the  entire  unstable  spectrum  of  the  current-driven  ion-cyclo¬ 
tron  modes,  it  is,  therefore,  necessary  to  solve  the  dispersion 
relations  numerically.  Indeed,  numerical  analysis  shows  that 
the  assumed  conditions  (|;c|  <C  1  and  :nj  'Jk>  1  for  all  n)  are 
rather  difficult  to  satisfy,  particularly  when  the  maximum 
values  of  the  growth  rates  occur. 


(ImW) 


(ImW) 


max 


FIG.  2.  ImW(=  Imco/Qj),  maximized  with  respect  to  k  and  k  .  as  a  func¬ 
tion  of  V0(=  V0/1/,||)  for  (ope/Sle  =  1/15  and  selected  values  of  7',, | /7)|| 
and  /,  7;  .  The  solid  curves  represent  results  for  bi-Maxwellian  plasma 
and  the  dashed  curves  represent  results  for  kappa-Maxwellian  ( ;c,.  —  tc,  —  3) 
plasma.  The  curves  are  labeled  by  the  selected  values  of  Te^/T^  and 
Tq/Tu. 


IV.  NUMERICAL  SOLUTIONS  OF  THE  DISPERSION 
RELATIONS 

We  have  developed  a  numerical  code  that  evaluates  the 
plasma  dispersion  functions  ZM(c),  ZKJc)  and  solves  the  dis¬ 
persion  relations  to  find  the  complex  values  of  co(=  cor  +  iy) 
as  a  function  of  the  aforementioned  variable  parameters. 
Here,  we  consider  the  electron-proton  plasma 
[me/mj=  1/1836)  and  concentrate  on  the  solutions  with 
CDr  «  Q,.  For  kappa-Maxwellian  plasma,  Ke  =  fc,  =  3  is 
assumed.  The  numerical  results,  graphically  presented  in 
Figs.  1-9,  13,  and  14,  are  obtained  by  assuming 

cope/ Qe  =  1/15,  which  is  consistent  with  the  assumption 
<))pe/Qr  -C  1  used  in  the  derivation  of  the  dispersion  rela¬ 
tions.  The  important  effects  of  cope/Sle  on  the  stability  prop¬ 
erties  are  discussed  later  in  this  section  by  means  of  Figs. 
10-12. 

Figure  1  shows  the  threshold  values  (critical  values)  of 
Vo  (normalized  to  Ve\\  and  denoted  by  Vc),  above  which  the 
instability  is  excited,  as  a  function  of  Te^/T^  and  T^/T^. 
The  threshold  values  are  obtained  by  minimizing  Vo/Vf|| 
with  respect  to  both  k  and  k ,  .  In  this  figure  and  in  all 
the  other  figures,  the  solid  curves  represent  the  results  for  bi- 
Maxwellian  plasma  and  the  dashed  curves  for  kappa- 
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(ReW) 

max 
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max 


FIG.  3.  ReW(=  Reco/Q/),  corresponding  to  (IrnW)max  shown  in  Fig.  2,  as 
a  function  of  V0(=  Vo/Ve||)  for  cope/Qe  =  1/15  and  selected  values  of 
T4/Til{  and  Tq/Tu_.  The  solid  curves  represent  results  for  bi-Maxwellian 
plasma  and  the  dashed  curves  represent  results  for  kappa-Maxwellian 
( Ke  —  Ki  =  3)  plasma.  The  curves  are  labeled  by  the  selected  values  of 

7Vr;||  and  Tin /r,j_. 


Maxwellian  plasma.  The  top  panel  of  Fig.  1  shows  that,  for  a 
fixed  value  of  Tj\\/Tj±,  the  value  of  Vc  in  both  types  of  plas¬ 
mas  decreases  sharply  as  Te^/T^  increases  and  then  settles 
down  to  a  value  which  is  nearly  independent  of  Tm/T^.  It 
also  indicates  that,  for  a  fixed  value  of  Te\\/Tj\\,  the  value  of 
Vc  decreases  as  Tm/T(±  increases.  This  is  shown  more  ex¬ 
plicitly  in  the  bottom  panel  of  Fig.  1.  The  results  for  bi- 
Maxwellian  plasma  are  in  agreement  with  those  obtained 
earlier.2 The  new  results  in  Fig.  1  are  that  a  larger  Vc 
(i.e.,  larger  current)  is  needed  to  excite  the  instability  in 
kappa-Maxwellian  plasma  for  all  Tm/Tu_  and  for  all  Te  |  /T, 
above  a  certain  value,  which  increases  with  decreasing 

Tq/Ti±. 

Figure  2  shows  the  maximum  values  of 
Im  W(=  Im  co/fl,-)  [maximized  with  respect  to  k\\,  k  and 
denoted  by  (Im  W)max  ]  as  a  function  of  V0(=  Vo/Ve||)  for 
selected  values  of  Te\\/T^  and  Tj\\/Tj±.  Figure  3  shows  the 
corresponding  values  of  Re  W(=  Re  rn/Q,),  denoted  by 
(Re  W)max,  as  a  function  of  V0(=  Vo/Fe||)  for  the  same  val¬ 
ues  of  Ten/Tj\\  and  Tm/Tu_.  As  expected,  for  all  values  of 
Tel\/Tq,  Tfn/Ta.  and  for  both  types  of  plasma,  (ImW)max 
decreases  as  V0  decreases  toward  the  threshold  values.  How¬ 
ever,  the  new  and  interesting  result  is  that,  for  all  values  of 
T^/T,  and  T,  \/TlL,  (ImW)max  in  kappa-Maxwellian  plasma 


is  larger  than  that  in  bi-Maxwellian  plasma  when  V0  is  much 
larger  than  the  threshold  value;  decreases  more  rapidly  as  V0 
decreases;  and  then  becomes  smaller  as  V0  approaches  the 
threshold  value  Vc-  This  is  consistent  with  the  results  pre¬ 
sented  in  Fig.  1,  namely,  Vc  is  larger  for  kappa-Maxwellian 
plasma.  Figure  3  shows  that  (Re  W)max  is  smaller  in  kappa- 
Maxwellian  plasma  than  in  bi-Maxwellian  plasma  for  all  val¬ 
ues  of  V0  >  Vc,  re||/r,-||,  and  7/  | /T,  L .  It  should  be  pointed 
out  that  the  approximate  analytical  solutions  in  Sec.  Ill  do 
not  show  the  dependence,  albeit  weak,  of  the  real  frequency 
on  V0.  As  Fig.  2  indicates  the  maximum  linear  growth  rate 
of  the  instability  can  be  somewhat  large  when  Te\\/T^, 
Ti»/Tj±  are  large  (>  1)  and  simultaneously  when  V0  is  also 
large  (e.g.,  V0  =  1).  This  is  particularly  true  for  the  kappa- 
Maxwellian  plasma.  The  linear  analysis  still  remains  valid  as 
(Im  W)max/ (Re  W)max  <  0.2  under  those  conditions.  All  it 
suggests  is  a  faster  onset  of  the  nonlinear  processes.  How¬ 
ever,  according  to  Fig.  1,  the  instability  can  be  excited  with  a 
much  smaller  V0(<C  1)  when  T^/T^  >  1  and  T^/T,±_  >  1. 
For  example,  (1)  when  Te^/T^  =  1  and  T^/Tj±  =  5, 
Vc,m  —  0.28  and  Vc.k  —  0.32;  (2)  when  Te\\/Tj\\  =5  and 
Ti\\/Ti±  =  1,  Vc,m  —  0.14  and  Vc,k  —  0.17;  and  (3)  when 
Te\\/Tj\\  =  5  and  T,t/Tu_  =  5,  VC,M 0.12  and  VC,K  =  0.14. 
Here,  Vc,m(k)  refers  to  the  threshold  value  Vc  for  bi- 
Maxwellian  (kappa-Maxwellian)  plasma.  The  maximum  lin¬ 
ear  growth  rates  of  the  instability  are  quite  small  when  V0  is 
above  and  near  these  small  threshold  values  (see  Fig.  2)  and, 


(ImW) 

max 


(Kratio) 

max 


Tepar/Tipar 
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FIG.  4.  (Top  panel)  ImW(=  Imco/fi;),  maximized  with  respect  to  iy  and 
k  i .  and  (bottom  panel)  ),  for  which  the  maximum  values  of  ImW  are 

obtained,  as  a  function  of  T^  /T^  for  cope/Cle  =  1/15,  V0(=  Vo/V^y)  =  1, 
and  selected  values  of  Ta/Ta_.  The  solid  curves  represent  results  for  bi- 
Maxwellian  plasma  and  the  dashed  curves  represent  results  for  kappa- 
Maxwellian  (tc,,  =  k,  —  3)  plasma.  The  curves  are  labeled  by  the  selected 
values  of  r,y  /Tn_ . 
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FIG.  5.  (Top  panel)  ImW(=  Imco/£2,),  maximized  with  respect  to  k«  and 
k  i ,  and  (bottom  panel)  (£|  jk  ),  for  which  the  maximum  values  of  ImW  are 
obtained,  as  a  function  of  r,-||/r,x  for  o)pe/Sle  =  1/15,  V0(=  Vo/V/,||)  =  1> 
and  selected  values  of  Te^/T^.  The  solid  curves  represent  results  for  bi- 
Maxwellian  plasma  and  the  dashed  curves  represent  results  for  kappa- 
Maxwellian  (  k€  =  K;  =  3)  plasma.  The  curves  are  labeled  by  the  selected 
values  of  re||/T;||. 


for  such  values  of  V0,  the  linear  stability  analysis  is  certainly 
appropriate.  In  the  following  presentation  of  the  various 
characteristics  of  the  unstable  ion-cyclotron  modes  in  the 
two  types  of  plasma,  we  have  used  V0  =  1  in  order  to  cover 
a  large  range  of  values  of  Te«/T^  and  Tm/T^  for  which  the 
instability  occurs  (see  Fig.  1). 

The  top  panel  of  Fig.  4  shows  (ImW)max,  defined  as 
above,  as  a  function  of  Te«  /T, -y  for  selected  values  of  Tm /7/y. 
It  is  evident  that,  for  a  fixed  7,  j  /T,  L ,  (ImW)max  increases 
with  re||/7’,1|  (rate  of  increase  being  larger  for  a  larger 
Ti\\/Ti±)  and  that  it  increases  more  rapidly  in  kappa- 
Maxwellian  plasma.  The  bottom  panel  shows  the  values  of 
k\\/k±  [denoted  by  (Kratio)max  ],  for  which  (ImW)max  are 
obtained,  as  a  function  of  7/ 1/7',  for  the  same  selected  val¬ 
ues  of  Tjti/Tix,  It  shows  that  the  angle  (between  k  and  Bo)  of 
propagation  of  the  maximum  unstable  modes  becomes 
smaller  with  increasing  values  of  Te^/T^  (i.e.,  the  modes 
propagate  more  parallel  to  B0)  and  that  this  change  in  the 
angle  of  propagation  occurs  more  rapidly  in  kappa- 
Maxwellian  plasma.  Figure  4  also  indicates  that,  for  a  fixed 
7’(,||/7’,||,  both  (ImW)max  and  (Kratio)max  increase  substan¬ 
tially  as  Tj«/Ta_  increases.  These  are  shown  more  explicitly 
in  the  top  and  bottom  panels  of  Fig.  5.  Numerical  analysis 
further  shows  that  the  magnitudes  of  (ImW)max  and 
(Kratio)max  as  a  function  of  Te»/T^  and  T^/Tn  in  both  types 
of  plasma  depend  on  the  choice  of  the  value  of  V0  >  Vc- 


Tepar7Tipar=.5,  1.  2.  5 


Tcpar/Tipar=.5,  I.  2.  5 
Kpar=.08;  Tipar/Tiperji=l 


FIG.  6.  (Top  panel)  ReW(=  Reco/Q,  )  vs.  K  i  (  =  Kperp)  —  k  \  pi  and  (bot¬ 
tom  panel)  ImW(=  Imco/Q,)  vs.  K  for  cope/Cle  =  1/15, 
V0(=  Vo/Te||)  =  1,  K||(=  Kpar)  =  k^pt  =  0.08,  Tq/Tj±  =  1,  and  selected 
values  of  Ten/T^.  The  solid  curves  represent  results  for  bi-Maxwellian 
plasma  and  the  dashed  curves  represent  results  for  kappa-Maxwellian 
(k(,  —  k j  =  3)  plasma.  The  curves  are  labeled  by  the  selected  values  of 
Te\\lT,\\- 


Figures  6-9  show  the  spectral  (m  vs  k)  behavior  of  the 
unstable  ion-cyclotron  modes  for  the  selected  values  of 
7/,||/7/||  and  T^/Tj±.  In  these  figures,  Kpar(=  K||)  =  k\\pt  and 
Kperp(=  K  )  =  k ;  ph  where  p,  is  the  ion  gyroradius.  Fig¬ 
ures  6  and  7  show  Re  W  and  ImW  versus  K  for  K  =  0.08 
and  various  values  of  Te^/T^  and  T/n/T/x.  They  show  that, 
for  a  given  T,  |  /T,  (or  T,,  |  /T,  \ ),  Re  in  moves  closer  to  Q,-  and 
the  growth  rate  (Imm)  decreases  (suggesting  larger  critical 
drift)  as  7’(,||/7’,  (or  T^/T,±)  decreases.  Conversely,  Rem 
moves  away  from  Q,  and  the  growth  rate  increases  (suggest¬ 
ing  smaller  critical  drift)  as  Tm/T^  (or  7/||/7/x)  increases. 
These  conclusions  hold  for  both  bi-Maxwellian  and  kappa- 
Maxwellian  plasmas.  However,  as  the  figures  show,  there  are 
significant  quantitative  differences  between  the  unstable 
spectra  in  the  two  types  of  plasma.  Magnitudes  of  these  dif¬ 
ferences  depend  on  the  values  of  Te  |  /T,  and  T,  |  /T,  for  the 
assumed  values  of  Ke,  k,,  and  V0.  The  real  frequencies  are 
smaller  in  kappa-Maxwellian  plasma  than  in  bi-Maxwellian 
plasma  for  all  values  of  K  .  On  the  other  hand,  the  growth 
rates  are  larger  in  kappa-Maxwellian  plasma  than  in  bi- 
Maxwellian  plasma  beyond  some  values  of  K  .  Addition¬ 
ally,  the  unstable  spectra  in  kappa-Maxwellian  plasma 
extend  to  comparatively  larger  values  of  K  .  Figures  8  and  9 
show  Re  W  and  ImW  versus  Ky,  for  K  i  =  0.8  and  values  of 
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FIG.  7.  (Top  panel)  ReW(=  Rero/£2,)  vs.  K  ,  (  =  Kpcrp)  =  k  \  p-  and  (bot¬ 
tom  panel)  ImW(=  Imco/O,)  vs.  K  for  a>pi,/Cle  =  1/15, 
V0(=  Vn/V',,||  )  =  1.  K|[(=  Kpar)  =  k  |  p,  =  0.08,  7).  |  /7)  =  1,  and  selected 
values  of  7",  /TL  .  The  solid  curves  represent  results  for  bi-Maxwellian 
plasma  and  the  dashed  curves  represent  results  for  kappa-Maxwellian 
( k,.  —  k,  =  3)  plasma.  The  curves  are  labeled  by  the  selected  values  of 

T||Ao, 


FIG.  8.  (Top  panel)  ReW(=  Reco/Q,)  vs.  K||(=  Kpar)  —  k\\pt  and  (bottom 
panel)  ImW(=  Imoj/Q,)  vs.  Ky  for  cope/Qe=  1/15,  V0(=  Vo/V^u)  =  1, 
Kx(=  Kperp)  =  k±p(  =  0.8,  T^/Tj±  =  1,  and  selected  values  of  Te^/T^. 
The  solid  curves  represent  results  for  bi-Maxwellian  plasma  and  the  dashed 
curves  represent  results  for  kappa-Maxwellian  (Ke  =  K\  =  3)  plasma.  The 
curves  are  labeled  by  the  selected  values  of  Te\\/Tj\\. 


Te\\/Ti\\  and  T^/Ti±  same  as  those  in  Figs.  6  and  7.  They 
show  similar  (to  Figs.  6  and  7)  functional  dependence  of  the 
unstable  spectra  on  Te«/T^,  and  similar  differences 

between  the  unstable  spectra  in  the  two  types  of  plasma  as  a 
function  of  Ky.  Both  the  K  i  -  and  the  K  -dependence  of 
Im  W  and  their  differences  in  the  two  types  of  plasmas  may 
be  understood  in  terms  of  the  combined  effects  of  the  elec¬ 
tron  drive  (Landau  growth)  and  the  ion  cyclotron  damping. 

For  heavier-ion  plasma,  qualitatively  similar,  but  quanti¬ 
tatively  different  due  to  increased  ion  mass,  results  are 
obtained.  We  do  not  present  those  results  here.  Instead,  more 
interesting  effects  of  the  parameter  cope/Q.e  on  the  stability 
properties  of  the  ion-cyclotron  modes  are  presented  next. 
Figure  10  shows  the  variation  of  Vc  (threshold  value  of  Vo 
normalized  to  Ve\\)  with  respect  to  ojpe/Cle  =  a  (Sigma)  for 
Tel\/Tq  =  1  and  different  values  of  T^/Ti±.  When  a  is  less 
than  0.02,  Vc  is  large  exceeding  unity.  As  a  increases,  Vc 
decreases  sharply  and  then  settles  down  to  a  value,  which  is 
nearly  independent  of  a.  Furthermore,  the  values  of  Vc  are 
larger  for  smaller  values  of  7’,  |  / T,  and  for  all  values  of  er. 
The  same  functional  behavior  is  found  for  both  types  of 
plasma.  However,  as  evident  from  the  figure,  there  are  appre¬ 
ciable  quantitative  differences  between  the  two  types  of 
plasma.  The  top  and  bottom  panels  of  Fig.  11,  which  are 


obtained  for  V0  =  1,  T^/T^  =  1,  and  selected  values  of 
Tj\\/Ti±,  show  that  (ImW)max  and  (Re  W)max,  defined  as 
above,  decrease  as  a  decreases.  The  rate  of  decrement 
increases  as  T^/Ti±  increases,  and  the  instability  disappears 
when  a  is  smaller  than  a  certain  value  that  depends  on 
Tj\\/Tjj_.  The  disappearance  of  the  instability  has  to  do  with 
the  fact  that  the  electron  drive  decreases  while  the  ion 
damping  increases  as  a  decreases.  Finally,  Fig.  12  shows  the 
dependence  of  (Kratio)max,  defined  as  above,  on  a  for 
V0  =  1,  T4/Ti{l  =  1,  and  selected  values  of  7)  |  /  T,- , .  For  both 
types  of  plasma,  the  angle  (between  k  and  Bo)  of  propagation 
of  the  maximum  unstable  mode  decreases,  i.e.,  the  mode 
propagates  more  parallel  to  Bo,  with  increasing  a.  Quantita¬ 
tive  differences  between  the  two  types  of  plasma  are  evident 
from  the  figures.  Similar  dependences  on  a  are  found  for  other 
values  of  T^/T^,  but  they  are  not  shown  here.  We  have  cho¬ 
sen  a  =  1  / 15  in  the  presentation  for  illustration  purposes  only 
as  Vc  tends  to  be  insensitive  to  a  around  this  value. 

V.  QUASILINEAR  ION  HEATING  RATES 

It  has  been  argued  that  quasilinear  plateau  formation  can¬ 
not  stabilize  the  current-driven  ion-cyclotron  instability 
because  of  the  slowly  convecting  nature  of  the  modes  and 
growth  continues  until  other  nonlinear  mechanisms  provide 
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FIG.  9.  (Top  panel)  ReW(=  Recn/Q,-)  vs.  K||(=  Kpar)  —  k\\pi  and  (bottom 
panel)  ImW(=  Imcu/Q,-)  vs.  Ky  for  (Dpe/Cle  —  1/15,  V0(=  Vo/^||)  —  1, 
Kx(=  Kperp)  =  kj_p,  =  0.8,  Te\\/Ti\\  =  1,  and  selected  values  of  Ti\\/Ti±. 
The  solid  curves  represent  results  for  bi-Maxwellian  plasma  and  the  dashed 
curves  represent  results  for  kappa-Maxwellian  (k€  =  Ki  =  3)  plasma.  The 
curves  are  labeled  by  the  selected  values  of  T/ii/T/l. 


22  23  33  •  r 

saturation.  '  Quasilinear  analysis,  therefore,  does  not 
yield  correct  estimates  of  the  ion  heating  rates.  Nevertheless, 
for  the  purposes  of  comparison,  we  present  here  the  quasilin¬ 
ear  estimates  of  the  ion  heating  rates  due  to  ion-cyclotron  tur¬ 
bulence  in  bi-Maxwellian  and  kappa-Maxwellian  plasmas. 

The  quasilinear  response  of  the  magnetized  ions  to 
unstable  electrostatic  held  fluctuations  is  described  by34 
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(26) 


for  small  Imcitk  >  0.  Here  /t  =  k±v±/Qi,  =  (o(k±,  h\ )  is 
the  complex  solution  of  the  linear  dispersion  relation,  and 
<f>( k,  t)  is  the  potential  huctuation,  which,  within  the  frame¬ 
work  of  the  quasilinear  theory,  evolves  adiabatically  in  time 
according  to 


||<KM)|2  =  2Im  tOk(0l<KM)|2-  (27) 


v 

C 


FIG.  10.  Threshold  (critical)  value  of  electron  drift  speed,  normalized  to  Ven 
and  denoted  by  Vc,  as  a  function  of  cope/Qe  —  o  (Sigma)  for  Te\\/Tj\\  =  1 
and  selected  values  of  Tjn/Tu_,  The  solid  curves  represent  results  for  bi- 
Maxwellian  plasma  and  the  dashed  curves  represent  results  for  kappa- 
Maxwellian  [k€  =  Ki  =  3)  plasma.  The  curves  are  labeled  by  the  selected 
values  of  Tq/Tu_. 


We  assume  that  the  fundamental  harmonic  is  the  dominant 
mode  and  so  keep  n—  1  term  only.  Then,  in  the  resonant 
region  of  velocity  space  (Reco^  —  k\\V\\  —  Q ,  =  0),  Eq.  (26) 
can  be  approximated  by34 


(TmW) 

max 


(ReW) 

max 


FIG.  11.  (Top  panel)  ImW(=  Imcu/Q,-),  maximized  with  respect  to  k\\  and 
k±  and  (bottom  panel)  ReW(=  Recu/Q,-),  corresponding  to  the  (ImW)max,  as 
a  function  of  cope/Qe  =  a  (Sigma)  for  V0(=  Vo/V^u)  =  1,  Te^/T^  =  1,  and 
selected  values  of  T^/Ti±.  The  solid  curves  represent  results  for  bi- 
Maxwellian  plasma  and  the  dashed  curves  represent  results  for  kappa- 
Maxwellian  {kc  =  Ki  =  3)  plasma.  The  curves  are  labeled  by  the  selected 
values  of  7/||/7Yj_. 
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(Kralio) 


Sigma 


FIG.  12.  Values  of  (k\\/k±),  for  which  the  maximum  values  of  ImW  shown 
in  Fig.  11  are  obtained,  as  a  function  of  a)pe/Q.e  =  c  (Sigma)  for 
V0(=  Vo/Ve||)  =  1,  T4/Ti{\  =  1,  and  selected  values  of  T^/Tj±.  The  solid 
curves  represent  results  for  bi-Maxwellian  plasma  and  the  dashed  curves 
represent  results  for  kappa-Maxwellian  ( Ke  —  Kj  —  3)  plasma.  The  curves 
are  labeled  by  the  selected  values  of  T4/Tj±. 
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under  the  small  Im  cok  assumption.  Here,  <)(x)  is  the  Dirac  S 
function.  Multiplying  Eq.  (28)  by  m,Vj  /2,  using  Eqs.  (2)  and 
(7)  for  Fi  on  the  right-hand  side  (for  estimation  purposes), 
and  integrating  over  velocities,  we  find 
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for  kappa-Maxwellian  plasma,  provided  k,-  >1/2.  Here 
>7,k  =  (Remk  —  Q/)/(\/2k||V/||).  Figure  13  shows  I(kx,£||)  as 
a  function  of  the  unstable  values  of  Kq (=  k||p;),  for 

Kx(=  k±Pj)  =  0.8,  cope/ a  =1/15,  VO(=Vo/Ve||)  =  l, 
r,||/2",||  =  1,  and  selected  values  of  7'1-||/7’l-x.  The  solid  curves 
represent  Ibm  given  by  Eq.  (30)  for  bi-Maxwellian  plasma 
and  the  dashed  curves  represent  Ikm  given  by  Eq.  (31)  for 
kappa-Maxwellian  plasma.  The  figure  shows  that  Ikm  >  Ibm 


FIG.  13.  The  function  I  appearing  in  Eq.  (29)  vs.  unstable  values  of 
K||(H*||pf)  for  K±(=  k±pl)  =  0.8,  cope/Slt=  1/15,  V0(=;  Vo/V^,,)  =  1, 
7',,  1 1  /l)  =  1,  and  selected  values  of  r,|i/r,x.  The  solid  curves  represent  I 
given  by  Eq.  (30)  for  bi-Maxwellian  plasma  and  the  dashed  curves  represent 
I  given  by  Eq.  (31)  for  kappa-Maxwellian  ( k,.  —  Kj  =  3)  plasma. 


over  almost  the  entire  unstable  Km  spectrum.  This  is  basically 
due  to  the  power-law  dependence  of  Ikm  on  f7,k’  in  contrast 
with  the  exponential  dependence  of  IBM  on  ;/,k,  when 
\rjik\  >  1.  Similar  results  are  found  for  other  values  of 
Te\\/Ti\\. 

Next,  multiplying  Eq.  (28)  by  m,-v|,  using  Eqs.  (2)  and 
(7)  for  Fj  on  the  right-hand  side  (for  estimation  purpose), 
and  integrating  over  velocities,  we  find 


d_  _  Imrij  e2Q, 

dt  111  _  V  T4 


where 


7r~r  |0(k,  t)|2Ai(fei)J(kx,k||),  (32) 


(33) 


J(*x,*||)  —  (Reojk/Q,  —  1)IBm  =  Jbm, 

for  bi-Maxwellian  plasma,  and 

J(*x,*||)  =  (Rett»k/I2,  —  1)Ikm  =  Jkm,  (34) 

for  kappa-Maxwellian  plasma.  Figure  14  shows  J(kx,k||)  as 
a  function  of  the  unstable  values  of  K||  (=  k«  pt),  for 
K±(=  k±Pi)  =  0.8,  Wpe/Fle  =  1/15,  VO(=Vo/Ve!|)  =  l, 
T eii /Tm  =  1,  and  selected  values  of  Thi/T/x.  The  solid  curves 
represent  JBm  given  by  Eq.  (33)  for  bi-Maxwellian  plasma 
and  the  dashed  curves  represent  Jkm  given  by  Eq.  (34)  for 
kappa-Maxwellian  plasma.  The  figure  shows  that  Jkm  >  Jbm 
over  almost  the  entire  unstable  Ky  spectrum.  Once  again, 
this  is  basically  due  to  the  power-law  dependence  of  Jkm  on 
i'll k,  in  contrast  with  the  exponential  dependence  of  JBm  on 
)7,k,  when  tjik \  >  1.  Similar  results  are  found  for  other  values 
of  Tf||/T,-||. 

Furthermore,  we  may  estimate  that 
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for  bi-Maxwellian  plasma  and 
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FIG.  14.  The  function  J  appearing  in  Eq.  (32)  vs.  unstable  values  of 
K||(=  fc||pf)  for  Kx(=  k±Pj)  =  0.8,  wpe/Cle=  1/15,  V0(=  V0/V,||)  =  1, 
7’,,||/7);  =  1,  and  selected  values  of  7)||  /7',  .  The  solid  curves  represent  J 

given  by  Eq.  (33)  for  bi-Maxwellian  plasma  and  the  dashed  curves  represent 
J  given  by  Eq.  (34)  for  kappa-Maxwellian  ( Ke  —  k,  =  3)  plasma. 

for  kappa-Maxwellian  plasma.  Although  the  right-hand  side 
of  Eqs.  (35)  and  (36)  is  formally  the  same,  they  are  quite  dif¬ 
ferent  for  bi-Maxwellian  and  kappa-Maxwellian  plasmas.  In 
fact,  numerical  solutions  of  the  linear  dispersion  relations 
show  that  (Rea>](  —  Q,)/Q,  <C  1  for  both  types  of  plasma, 
but  (Re  04  —  Q,-) /Q,  has  smaller  values  for  kappa- 
Maxwellian  plasma,  when  Vo,  Tm/T^,  and  T\\/Tj±  are  not 
too  large  [see  Figs.  6-9].  Equations  (35)  and  (36)  then  sug¬ 
gest  that,  under  these  conditions,  ions  are  heated  preferen¬ 
tially  in  the  perpendicular  direction  and  this  effect  is 
enhanced  in  kappa-Maxwellian  plasma.  These  conclusions 
can  also  be  derived  from  Figs.  13  and  14  where  IBM,  IKM, 
JBM,  and  Jkm  are  shown. 

VI.  SUMMARY 

We  have  presented  for  the  first  time  the  linear  stability 
properties  of  the  current-driven  electrostatic  ion-cyclotron 
modes  in  kappa-Maxwellian  plasma  where  the  parallel  ve¬ 
locity  distributions  of  electrons  and  protons  are  modeled  by 
kappa  distributions  and  the  perpendicular  velocity  distribu¬ 
tions  are  modeled  by  Maxwellian  distributions  with  parallel 
temperatures  different  from  perpendicular  temperatures.  In 
particular,  dependence  of  the  stability  properties  on  various 
parameters,  such  as  Vo,  k\\,  k±,  Te«/Tj«,  Tm/Tu_,  and  <Dpe/kle , 
has  been  discussed  in  some  details.  For  comparison,  we  have 
also  presented  the  corresponding  results  for  bi-Maxwellian 
plasma  where  both  the  parallel  and  the  perpendicular  veloc¬ 
ity  distributions  are  modeled  by  Maxwellian  distributions 
with  parallel  temperatures  different  from  perpendicular  tem¬ 
peratures.  Quasilinear  estimates  of  resonant  ion  heating  rates 
due  to  ion-cyclotron  turbulence  in  the  two  types  of  plasma 
have  also  been  presented  for  comparison.  In  order  to  gain 
some  preliminary  understanding  of  the  linear  instability  in 
kappa-Maxwellian  plasma  and  its  comparison  with  that  in 
bi-Maxwellian  plasma,  we  have  included  approximate  ana¬ 
lytical  solutions  of  the  dispersion  relations.  But,  more  em¬ 
phasis  has  been  given  on  the  numerical  solutions  of  the 
dispersion  relations  in  order  to  provide  accurate  and  compre¬ 
hensive  comparison  over  a  wide  range  of  values  of  the  vari¬ 
ous  plasma  parameters.  Numerical  analysis  has  been  carried 


out  for  electron-proton  plasma.  Qualitatively  similar,  but 
quantitatively  different  due  to  increased  ion  mass,  results  are 
expected  for  heavier-ion  plasma. 

Linear  stability  properties  of  the  current-driven  electro¬ 
static  ion  cyclotron  modes,  as  well  as  the  quasilinear  ion 
heating  rates  due  to  ion-cyclotron  turbulence,  in  bi- 
Maxwellian  and  kappa-Maxwellian  plasmas  are  qualitatively 
similar.  However,  there  are  quantitative  differences,  which 
can  be  significant,  depending  on  the  values  of  Ke  and  ;c,.  The 
important  conclusions  of  the  comparative  study  are  summar¬ 
ized  as  follows. 

(1)  The  threshold  (critical)  value  of  the  electron  drift  speed 
(i.e.,  current),  above  which  the  current-driven  ion-cyclo¬ 
tron  instability  is  excited,  is  larger  in  kappa-Maxwellian 
plasma  for  all  values  of  /T,±  and  for  all  Te«/T^  above 
certain  value,  which  increases  with  decreasing  /Ttl 
(see  Fig.  1). 

(2)  Above  and  near  the  threshold  (critical)  value  of  the  elec¬ 
tron  drift  speed,  the  maximum  growth  rates  are  smaller 
in  kappa-Maxwellian  plasma.  But,  as  Vo  is  increased,  the 
maximum  growth  rates  become  larger  in  kappa- 
Maxwellian  plasma.  Additionally,  the  unstable  spectra  in 
kappa-Maxwellian  plasma  extend  to  comparatively 
larger  values  of  Ky(=£x/9,)  and  Km(=  k\\pt)  (see 
Figs.  2,  6-9).  The  frequencies  of  the  unstable  modes  are 
smaller  in  kappa-Maxwellian  plasma  (see  Figs.  3,  6-9). 
These  conclusions  are  valid  for  all  values  of  Te\\/Tj\\  and 
Ti\\/Ti±  for  which  the  instability  can  be  excited  with  a 
given  Vo. 

(3)  The  angle  (between  k  and  Bo)  of  propagation  of  the  max¬ 
imum  unstable  mode  decreases  with  increasing  T^/T^ 
and  Tj\\/Tj±  in  both  bi-Maxwellian  and  kappa- 
Maxwellian  plasmas.  In  other  words,  as  Te»/T^  and 
Ti\\/Tj±  increase,  the  maximum  unstable  mode  propa¬ 
gates  more  parallel  to  Bo.  However,  the  change  in  the 
angle  of  propagation  occurs  more  rapidly  in  kappa- 
Maxwellian  plasma  (see  bottom  panels  of  Figs.  4  and  5). 
The  magnitude  of  the  change  in  the  angle  of  propagation 
as  well  as  the  magnitude  of  the  differences  between  the 
bi-Maxwellian  and  kappa-Maxwellian  plasmas  depends 
on  the  value  of  V0  >  Vc. 

(4)  The  threshold  (critical)  value  of  the  electron  drift  speed, 
above  which  the  current-driven  ion-cyclotron  instability 
is  excited,  decreases  sharply  as  (ope/£le  increases  and 
then  settles  down  to  a  value,  which  is  nearly  independent 
of  cope  /  Qr ,  for  both  bi-Maxwellian  and  kappa- 
Maxwellian  plasmas  (see  Fig.  10).  For  both  types  of 
plasma,  the  maximum  value  of  Im  W(=  Imco/Q,)  and 
the  corresponding  value  of  Re  W(=  Reco/Q,)  decrease 
as  cope/Q.e  decreases.  The  rate  of  decrement  increases  as 
T;  |  /Ti  increases,  and  the  instability  disappears  when 
cope/kle  is  smaller  than  a  certain  value,  which  depends  on 
Ti\\/Ti±  (see  Fig.  11).  For  both  types  of  plasma,  the  angle 
(between  k  and  Bo)  of  propagation  of  the  maximum 
unstable  mode  decreases,  i.e.,  the  mode  propagates  more 
parallel  to  Bo,  with  increasing  (ope/Qe  (see  Fig.  12). 
Quantitative  differences  between  the  two  types  of  plasma 
are  evident  from  the  figures. 
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(5)  Quasilinear  estimates  of  the  resonant  ion  heating  rates  sug¬ 
gest  that  both  the  perpendicular  and  the  parallel  heating 
rates  are  larger  in  kappa-Maxwellian  plasma.  This  is  basi¬ 
cally  due  to  the  power-law  dependence  (as  opposed  to  the 
exponential  dependence)  on  vy  of  the  distribution  of  reso¬ 
nant  ions  in  kappa-Maxwellian  plasma  [see  Eqs. 
(29) — (34)].  They  also  suggest  that,  under  certain  conditions, 
ions  are  heated  preferentially  in  the  perpendicular  direction 
in  both  bi-Maxwellian  and  kappa-Maxwellian  plasmas,  but 
this  effect  is  enhanced  in  kappa-Maxwellian  plasma. 

The  present  study  is  based  on  nonrelativistic  kappa  dis¬ 
tributions.  For  some  space  plasma,  highly  energetic  particles 
are  better  modeled  by  relativistic  kappa  distributions.  The 
reader  is  referred  to  the  works  by  Xiao  and  Xiao  et  al: 
for  relativistic  kappa  distributions. 
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